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Abstract 

Using the closed-time path integral approach, we nonperturbatively study inelastic tunneling of 
electrons via magnetic impurities in the barrier accompanied by phonon emission in a magnetic 
tunnel junction. The spectrum density of phonon emission is found to show a power-law infrared 
singularity ~ w"^^"^-' with g the dimensionless electron-phonon coupling. As a consequence, the 
tunneling conductance GiV) increases with bias voltage \V\ as GiV) — G(0) ~ |V^P^, exhibiting a 
discontinuity in slope at ^ = for g < 0.5. This theory can reproduce both cusp-like and non- 
cusp-like feature of the zero-bias anomaly of tunneling resistance and magnetoresistance widely 
observed in experiments. 
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The tunneling magnetoresistance (TMR) effect in a magnetic junction composed of two 
ferromagnetic (FM) electrodes separated by an insulating layer is presently of great interest 
because of its potential of device applications such as random access memories and magnetic 
sensor o'^i^i^i^i^i^i'^'i^i^i^^i^'^ . While the TMR effect has been attributed to spin polarization of 
conduction electrons in the FM electrodes^*^, its novel voltage dependence is still puzzling. 
The tunneling resistance and TMR ratio generally decrease with increasing bias voltage, 
often exhibiting a cusp-like feature at zero bia o^i^i^i'^'i^i^i^^i^'^ . This effect, usually called the 
zero-bias anomaly (ZBA)^, was not only widely observed in transition-metal magnetic junc- 
tions, but also exists in semiconductor GaMnAs/AlAs/GaMnAs magnetic junctionsi^. 

It is found that the ZBA is very sensitive to the insulating material of the barrier. In 
earlier measurements^ on Fe/Ge/Co junctions, the characteristic voltage V1/2 that halves 
the TMR ratio is about 3 mV. It was reported recently that for FM/AI2O3/FM with Fe, 
Co, Ni, or their alloys as the FM electrodes, V1/2 is usually one hundred to several hundred 
mV'^'^'2^, and for Ni/NiO/Co junctions V1/2 is a few ten mV— . Junctions made of the same 
FM electrodes and identical but differently prepared insulator vary considerably concerning 
the voltage dependence of the resistance and TMR—. It was reported that the ZBA in 
the semiconductor magnetic tunnel junctions is much more striking with V1/2 about 3m V—. 
Interestingly, in a very recent experimentiS, of the TMR between a Co sample and a CoFeSiB 
tip through a vacuum barrier, the ZBA did not occur for bias voltage upto a volt. 

Understanding of the ZBA is not only important for optimizing the TMR in practical 
applications, but also provides critical insight into the underlying physics of spin- dependent 
transport. There exist several possible explanations for the anomaly. The simplest one may 
involve the energy dependence of the electron density of states (DOS) and elastic trans- 
mission matrix elements^^i^^. Unfortunately, by comparison with experimental data, the 
TMR ratio calculated from such a model decreases too slowly at low bias^. The second 
mechanism is that the hot electrons from the emitting electrode may be scattered by lo- 
cal magnetic moments at the interfaces through s-d interaction^, which contributes to the 
conductance a term AG ~ \V\ at low bias voltage V. This theory succeeded in explaining 
some experimental data. The third mechanism might be based upon the electron tunnel- 
ing assisted by magnons and phonon o'^i^^i^'^i^^ , which was found to contribute AG ~ |V^p/^ 
for FM magnons, for surface antiferromagnetic magnons, and l^l"^ for phonons. This 
mechanism may explain some experimental data without the cusp-like feature at zero bias^^, 
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but has the difficulty to interpret the cusp-hke feature of the ZBA observed in many other 
experimenta^i^'^'iS^. On contrary to these intrinsic mechanisms, the effect of extrinsic non- 
magnetic and magnetic impurities was emphasized in many works^ i^^i^-^i^^i-^^i-^^i^^i^^i^^i^^ . The 
vohage dependence caused by the elastic electron tunneling via impurities was studied the- 
oretically in a few works^iJ^^. However, none of them reproduced the essential features of 
the ZBA. Based upon the results of an earlier theory-, some author o^°i^^ argued that, in 
the presence of inelastic scattering, while single-impurity tunneling processes are still elas- 
tic, electron tunneling via multiple-impurity chains may contribute an inelastic power-law 
conductance AG ~ \V\'^ with p = 4/3 (5/2) for two-impurity (three-impurity) processes. 
Presently, the relevant mechanism for the ZBA is still under debate. 

In this work, a combined effect of the magnetic impurities in the barrier and the electron- 
phonon interaction is proposed as a possible mechanism for the ZBA. What we consider in 
this model is not a simple addition of the elastic tunneling of electrons via the impurities 
and the inelastic one assisted by phononsi^, but a novel inelastic electron tunneling due 
to the strong interplay between them, in which the electron tunneling via impurities is 
always accompanied by phonon emission. It is found that the spectrum density of phonon 
emission has a power-law infrared singularity ~ lu^^^^^^ with g the dimensionless electron- 
phonon coupling. As a consequence, the conductance via impurities ^(V^) obeys a power- 
law voltage dependence G^{V) ~ with a discontinuity in slope at ^ = for g < 0.5. 
The conductance exponent 2g may be either smaller or greater than unity, being different 
from the earlier theorie o^'^^'^'^i^^i^^ where the conductance exponents are never smaller than 
unity. The calculated results can reproduce both cusp-like and non-cusp-like features of 
ZBA observed in experiments. An experiment based upon the isotope effect is suggested to 
further verify this theory. 

The model Hamiltonian of a magnetic tunnel junction is written as H = + H^^ + Hq + 
H^. Here, Hi^(r) = J^ks^kT^fhT^'^^kT^ the Hamiltonian of the electrons in the left (right) 
electrodes. The third term is the Hamiltonian of the electrons on the magnetic impurities. 

Ha = E(^o<5ss' - J(Tss' ■ Si)dldis> + ^ ujqxblxbqx + E ^.Ae^'^ '^' (&,a + &1,a) djsdis • (1) 

Iss' qX qXls 

Here, dis annihilates an electron with spin s on magnetic impurity /, and hq\ annihilates an 
acoustic phonon of wave vector q and polarization A in the system. The magnetic impurities 
in the barrier are presumably the atoms of FM metals of the electrodes mixed into the 
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barrier. Some of their outer-shell electrons are itinerant, and the others form localized 
spins. J represents the s-d like interaction between the itinerant electrons and localized 
spin Si with a the Pauli matrix. The spins of the transition-metal ferromagnets under 
consideration are usually greater than 1/2, and can be treated as classical approximately. 
Since the impurities are also a part of the lattice, electrons scattering with the impurities 
will generally cause emission or absorption of phonons of the system, which is described 
by the electron-phonon interaction of the last term in Eq. (0). It is worth mentioning that 
a similar impurity-mediated electron-phonon interaction has been used to account for the 
leading correction to residue resistivity of simple metalst^^^. Here, the direct scattering 
of electrons by phonons is neglected, which was found to open new tunneling channels, and 
yield even, positive correction to the conductance^^. It does not account for the ZBA. The 
electron tunneling Hamiltonian is given by 

= E {T^k'sO^M's + h-c.) + E {Tlld^ldi, + h.c.) + E {Tll4ldis + h-c) , (2) 

kk's Iks Iks 

where the first term represents direct tunneling of electrons between the electrodes, and 
the second and third terms represent transfer of electrons between the electrodes and the 
impurities. 

To treat the electron-phonon interaction in a nonperturbative manner, the Lang-Firsov 
transformation e'^ with S = — '}2q\is{^q\/^q\)Asdise^^'^'-{^q\ — ^'^-q\Y^ is performed to diago- 
nalize the electron-phonon interaction. The transformed Hamiltonian e^He~^ is the same as 
the original one, except that the last term in Hq vanishes, and operators dis and d\^ in are 
replaced with disXi and d\^X} , respectively, with Xi = exp[J2qx{Mgx/ujgx)e"^'^'{bgx - b^-qx)]- 
In this nonperturbative treatment, we see clearly that electron tunneling via the impurities 
causes phonon emission or absorption, which reversely renormalizes the electron tunneling 
matrix elements via the impurities. 

The path-integral formulism on the closed-time patb^S*^ is employed to study this 
nonequilibrium problem. By integrating out the electron variables in the two electrodes, the 
retarded (advanced) Green's function for electrons on impurity / is obtained in frequency 
domain as Gf^ria)i^) — [i^ — Eq-\- Ja ■ Si — T,i^r{a){^)]^^ ■ Here, 2x2 matrix representation in 
electron spin space has been employed. The retarded self-energy f^i^ri^^) = ^]^,ri^) + ^f,ri^) 
is a diagonal matrix with diagonal elements given by 

= IpT'IT^' J duj'{[l - fun.iu^ + u;')]B^{-u') + U,,,{u; + u;')B^{u;')} , (3) 
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and the advanced self-energy S;a(ci;) = [S;^,.(cj)]'''. Here, as usual, the momentum and spin 
dependence in the transmission matrix element T/^j^^ is neglected^, p^'^' and /L(R)(a;) are 
the electron DOS and Fermi distribution function of the left (right) electrode, respectively. 
By (to) is the Fourier transform of — z(X;(t)X/(0)), whose general expression can be found 
in reference^^. 

The tunneling current / is calculated from the rate of the charges flowing out of the 
emitting electrode. The total conductance G(y) = I /V is obtained as G{V) = G° + G^(y), 
where G° = 27re^|T°pTr(p'"p^) is the direct tunneling conductance independent of bias 
voltage, and 

G\V) = JdeJ dco, I d^25>K)5>(^2){/L(e + ^i)[l-/R(e-^2)] 



- Ue + u,) [1 - A(e - u,)]}\Tl-Tnm [p-Gf^MP^'Gf^ai^)] , (4) 

comes from the tunneling via A^i impurities. The bias dependence of the tunneling con- 
ductance is embodied in G^{V). In Eq. Q, the last factor is averaged over the po- 
sitions and spin orientations of the impurities, p^ = [(p| + pj) + (p| — pi)d'z]/2, and 

— [(Pt + Pi) ='= (Pt ~ Pi)^2;]/2 with sign + for parallel alignment and — for antiparal- 
lel alignment. The two electrodes are considered of the same material with electron DOS p| 
and pi for majority and minority spins, respectively. 

There are two distinct resonant energies E_ = Eq — JS and = Eq + JS, at which 
the tunneling conductance given in Eq. (0)) will be substantially enhanced. In this work, 
we confine ourselves to the off-resonance case, which is more probable in reality. The en- 
ergy differences \Ep — E±\ with Ep the Fermi level are considered to be much greater than 
the broadenings of the resonant states. Besides, the bias voltage V is also relatively small 
\eV\ < \Ep — E±\. In this case, S;^r(a)(e) in G^«Jr(a)(^) "^^^ ^e neglected. The last factor in Eq. 
(jU is then divided into two parts: {Tl^T™]"^ and TY[p^Gf^{e)p^Gf^{e)]. The former is a con- 
stant after averaged over the impurity position R;, and the latter involves only the average 
over the orientation of Si. The energy of the magnetic field used to produce the parallel align- 
ment of the FM electrodes, usually several ten or hundred oersteds {gi^p^SB < O.OlmeV), 
is much smaller than the entropy gain due to randomizing an impurity spin even at exper- 
imentally low temperatures, e.g., IK [A;BTln(2S'+ 1) ~ O.lmeV]. Therefore, it is assumed 
that the impurity spins are randomly oriented for both parallel and antiparallel alignments. 
Moreover, one can verify that the results will not change significantly, if the magnetic field 
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FIG. 1: Normalized impurity conductance of nonmagnetic tunnel junctions as a function of bias 
voltage at zero temperature. Here J = and p| = Pi- 

is strong enough to align all the spins ferromagnetically in the parallel alignment. 

The function B^{uj) represents the spectrum density of phonon emission and absorption. 
For simplicity, three-dimensional Debye's model is used to describe the acoustic phonons. 
A dimensionless electron-phonon coupling constant g is defined as (7 = 3 X^^a \Mq\\^ujq\/uj'^ 
with uoj^ the Debye frequency. Thus, at zero temperature we have 



where Cin(a;) = dt[l — cos(t)]/t and Si(x) = Jq dtsin{t)/t are the cosine and sine integral 
functions^^. G\V) via impurities calculated by use of Eqs. (4) and (0) is plotted in Fig. 1 
for nonmagnetic tunnel junctions. Since \img_,Q[iB^{Lj)] = 2n6{uj), \img^QG\V) = Gq is a 
constant independent of V, which is actually the maximum of G^{V). From Fig. 1, it is clear 
that G^{V) vanishes at V = for any nonzero g, indicating that tunneling via impurities is 
totally suppressed, and G\V) increases as the voltage |^| is elevated. The stronger is the 
electron-phonon coupling g, the slower the conductance G^{V) increases with \V\. 

To get some insight into low bias behavior of the tunneling conductance, we derive the 




dtexp {iut — g [Cm^Uot) + iSi^Uut)]} 



(5) 
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low-cij approximation of Eq. (0) 

iuJuBy{u) ^ ae{uj){uj/ujj,)-^'^-^^ , (6) 

where a = 2r(l — g) sm{gTT)e''^'^ with 7 = 0.5772 the Euler constant, and T{z) = 
dtt^~^e~^ is the ordinary gamma function^^. The unit step function 6{uj) indicates 
that at zero temperature only emission of phonons is possible. Substituting Eq. (jH)) into 
Eq. (jH), we obtain an interesting result G\V) ~ For g < 0.5 {g > 0.5), its slope 

dG^{V)/dV ~ \ V\'^^/V is discontinuous (continuous) at ^ = 0. 

In Fig. 2, the resistances for parallel alignment (Rp) and antiparallel alignment (-Rap) 
as well as TMR = (-Rap — -Rp)/-Rap of magnetic tunnel junctions are shown as functions 
of the bias voltage. It is found that -Rap always decreases faster than -Rp, resulting in a 
decline of TMR. This behavior can be understood by the following argument. At ~ 0, the 
direct tunneling dominates. For antiparallel alignment, the majority-spin (minority-spin) 
band in the left electrode is only connected to the minority-spin (majority-spin) band in the 
right electrode, leading to a relatively high resistance. With increasing \V\, the impurity 
tunneling increases. Electron hopping via the magnetic impurities causes a mixing of the two 
spin channels, which effectively connects up the majority-spin bands in the two electrodes. 
As a consequence, the resistance drops rapidly. For g = 0.4, the TMR is discontinuous 
in slope at = 0. According to Fig. 2, either increasing the electron-phonon coupling g 
or decreasing relatively the tunneling via impurities weakens the ZBA. In this theory, the 
energy scale of the ZBA is from zero to some u^,. In view of that in most materials 
ranges from lOmeV to lOOmeV, the characteristic energies of the ZBA observed in tunnel 
junctions with Ge^, NiO^^ and AlAs^^ as insulating barriers are consistent with our theory. 
For FM/AI2O3/FM, Vi/2 is greater than 100me V^'^i^'7i?i? . The impurity density might be low 
in these later junctions so that the ZBA is not prominent. According to a simplified model 
for impurity tunneling^^, we estimate that the impurity effect and the ZBA will become 
observable as the impurity number per unit cross section reaches 1 / fim'^ or higher. 

The present theory as well as the earlier theories^'i^'i^'i^^ all predict that G{V) for the 
tunnel junctions has a power-law dependence on small bias voltage, i.e., G(y) —G{0) ~ \ V\^. 
The low-temperature curves of the conductance (resistance) and TMR exhibit clear cusp- 
like feature at zero bias in most experimentsi^'^'i'^^^, and non-cusp-like feature in a few 
experiments^^. Mathematically, the cusp-like (non-cusp-like) feature requires that the slope 
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FIG. 2: Calculated results for i?p, -Rap) and TMR as functions of bias voltage for different g. The 
resistances are normalized by 2Rj^p{y = 0). Here, £^0 = E^-, and spin polarization [p-^ — Pi)/{p^ + 
Pi) = 0.6. 

of the conductance dG(y)/dV ~ \V\p/V be discontinuous (continuous) at ^ = 0, and so 
p < 1 {p > 1). As mentioned in the introduction, all the earlier theoriesi^'i^'i^'^i predicted 
constant p > 1, except for the theory of Zhang et al^, where p = 1. In the present 
theory, p = 2g may be either smaller or greater than unity, depending on the strength of 
the impurity-mediated electron-phonon interaction, which is more flexible for understanding 
both the cusp-like and non-cusp-like features within a unified theoretical framework. Our 
theory might be further verified directly by measurements of isotope effect. Let us assume 
that the insulating barrier is made of certain oxide. With gradual replacement of the normal 
O^^ atoms with isotope O^^ atoms, the mass M in a lattice cell increases, and the phonon 
frequency and Debye frequency change as tOqx ~ M~^/^ and ~ M~^/^. We notice 



8 



l-^gAp ~ {MiLjgx)~^ ~ M^/^, where Mi is the impurity mass rather than the cell mas o^^'^^ . 
Consequently, the conductance exponent 2g in the present theory will increase continuously 
as 2g ~ M'^/^ with the oxygen isotope enrichment. In contrast, the constant conductance 
exponents predicted by the earlier theories are expected not to change in this case. 
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